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Abstract. We report on recent results from VISH2+1, a code that solves the 
relativistic Israel-Stewart equations for causal viscous hydrodynamics for heavy- 
ion collisions with longitudinal boost invariance. We find that even "minimal" 
shear viscosity r\l s = h/(4n) leads to a large reduction of elliptic flow compared to 
ideal fluid dynamics. We explore systematically the sensitivity of this reduction 
to the equation of state, system size, initial conditions, and the microscopic 
relaxation time in different formulations of the Israel-Stewart equations. 



1. Introduction 

A good motivation for studying relativistic hydrodynamics for viscous fluids in order 
to describe the expansion of the fireballs created in relativistic heavy-ion collisions 
is given in Chaudhuri's plenary talk pQ. We here describe results from systematic 
studies done with VISH2+1, our numerical code for solving viscous Israel-Stewart 
[5] hydrodynamics for longitudinally boost-invariant systems in two transverse space 
dimensions. We will concentrate on the most important results and refer the reader 
to the original publications J3J HI El E] for details. We include only the effects of shear 
viscosity, neglecting bulk viscosity and heat conduction. 

2. Israel-Stewart formalism 

The Israel-Stewart formalism [2] resolves problems of acausality and numerical 
instbility of the relativistic Navier-Stokes equation by allowing the shear viscous 
pressure components t:^ v of the energy momentum tensor to evolve dynamically, on 
a short microscopic time scale 7V, towards their Navier-Stokes limit. Different forms 
of the equations governing the evolution of t:^ u are being discussed in the literature 
(see references in [fl |Hj). Although they all converge to the same Navier-Stokes limit 
as — > 0, their results disagree for finite r w (which is necessary for stable numerical 
implementation) [U |H [6l [7] ■ Recent work [6] has shown that the weakest sensitivity 
of physical observables to the "regulator" tv is obtained by employing the so-called 
"full" Israel-Stewart equations used in [6j[7l[8] whereas "simplified" versions advocated 
in [U [3] that are obtained by dropping certain terms in the "full" equations and 
were used in [TJ [H [5J [S] give results that depend strongly on the (presently poorly 
constrained) value for t„. As shown in [3] at the end of section IV, the extra terms 
in the " full I-S" formulation cause a reduction of the effective kinetic relaxation time 
for tt^ , especially in rapidly expanding systems. This keeps the viscous pressure 
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smaller than in the "simplified I-S" formalism, holding it always close to its Navier- 
Stokes limii[§] and thereby eliminating large excursions of n' 11 ' that increase viscous 
entropy production and lead to a stronger viscous reduction of the elliptic flow as r T 
is increased [5J. The Ottinger-Grmela framework used in [10) appears to give very 
similar results as the "full I-S" one[jj [IT]. 



3. Cooling rates for ideal and viscous heavy-ion fireballs 

In the early expansion stage of a heavy-ion collision, shear viscosity leads to a reduction 
of the longitudinal and an increase of the transverse pressure. The reduced longitudinal 
pressure decreases the work done by longitudinal expansion, thereby reducing the 
initial cooling rate of the fireball. This leads to a somewhat increased lifetime of the 
quark-gluon plasma phase. The increased transverse pressure, on the other hand, 
causes larger transverse acceleration and stronger radial flow of the matter than in 
ideal fluid dynamics. Due to the larger transverse flow, the center of the viscous fireball 
cools more quickly during the late stages than an ideal fluid, thereby slightly reducing 
the total fireball lifetime until freeze-out for central and near-central collisions (see 
Fig. [1]). The viscous fireballs created in peripheral collisions don't live long enough 
for this mechanism to manifest itself; they live longer than their ideal counterparts, 
due to the decreased initial cooling rate arising from the smaller longitudinal pressure 
(left points in Fig. [J). 
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Figure 1. (Color online) Central freeze-out times for Au+Au and Cu+Cu 
collisions with differential initial peak energy densities eo = e(0, 0; 6=0) as 
a function of final charged hadron multiplicity density. Points belonging to 
one series correspond to different impact parameters, with central (peripheral) 
collisions corresponding to large (small) multiplicity densities. Left: ideal fluid 
dynamics. Right: viscous fluid dynamics with rj/s = 0.08 and tv = 3r]/sT, using 
the full I-S equations. 



For an equation of state (EOS) with a first-order quark-hadron phase transition, 
the mixed phase matter is free of pressure gradients and hence not accelerated. 
Velocity gradients in the mixed phase generate, however, viscous pressure components 

§ We thank D. Rischke and U. Wiedemann for an illuminating discussion on this point. 
|| It has recently been checked that VISH2+1 with "full I-S" equations and the codes developed by 
P. & U. Romatschke [7] and Dusling & Teaney 1101 give, for identical initial conditions and equations 
of state, the same time evolution for the momentum anisotropy in non-central Au+Au collisions at 
RHIC energies (up to small numerical errors). We thank K. Dusling for providing us with his results 
for comparison, and P. Romatschke for help with using their code. 
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whose gradients continue to accelerate the viscous fluid even while it passes through 
the mixed phase where the thermal pressure gradients vanish. As a result, the viscous 
fluid spends less time in the mixed phase than the ideal one. 

In ideal hydrodynamics, a first-order phase transition with a mixed phase (MP) 
generates large velocity gradients near the QGP-MP and MP-HG (HG = hadron 
gas) interfaces. In viscous fluids such gradients generate viscous pressures that act 
against building up large velocity gradients. In consequence, all prominent structures 
in ideal hydrodynamics that arise from discontinuities in the speed of sound in ideal 
hydrodynamics are washed out by viscous effects. Shear viscosity thus effectively turns 
a first-order phase transition into a smooth cross-over transition. 

4. Momentum anisotropy in non-central collisions 

IrpXX rpy y \ 

In viscous hydrodynamics the total momentum anisotropy e p = na s two 

contributions. The first arises (as in ideal fluids) from anisotropies in the collective 
flow velocity. The second is due to anisotropies of the viscous shear pressure tensor 
which contribute a momentum anisotropy even in the local rest frame of the fluid. 
In heavy-ion collisions, these two contributions are found to be of opposite signs: the 
viscous pressure effects reduce the total omentum anisotropy. In the early stages 
of the expansion, the collective flow anisotropies evolve very similarly in viscous 
and ideal fluids; however, due to the large longitudinal expansion rate ~ 1/r, the 
viscous pressure effects are large at early times, leading to a strong reduction of 
the total momentum anisotropy. At late times, the viscous pressure effects are 
dramatically reduced, but by then their time- integrated effect has led to a slower build- 
up of collective flow anisotropies compared to ideal fluids. So the total momentum 
anisotropy at late times is again much smaller for the viscous fluid, even though the 
negative viscous pressure contributions to this anisotropy are now almost neglegible. 
Overall, shear viscous effects are seen to cause during all evolution stages a significant 
reduction of the momentum anisotropy compared to ideal fluids. 

5. Spectra and elliptic flow 

The larger radial flow caused by the viscous increase of the transverse pressure 
gradients reflects itself through flatter transverse momentum spectra for the finally 
emitted hadrons. To obtain similar spectral slopes in viscous and ideal hydrodynamics 
thus requires a retuning of initial conditions and/or freeze-out conditions. Starting the 
viscous expansion later with lower initial energy densities works [7J , but this neglects 
transverse flow generated during the earlier pre-equilibrium stages. A consistent 
retiming of initial conditions that is compatible with dynamical models for the pre- 
hydrodynamic stage remains an open issue. 

The smaller momentum anisotropies of viscous fireballs discussed above manifest 
themselves in a reduction of the elliptic flow coefficients extracted from the final hadron 
spectra. Initial elliptic flow results from our group [H [5] and Romatschke [7J seemed 
to show large disagreements: where we saw for b = 7fm Cu+Cu collisions at pr = 
2 GeV/c a reduction of pion V2 by about 63% from the ideal fluid value, Romatschke 
saw for b = 7fm Au+Au collisions only a 25% reduction. These differences have been 
understood and resolved in [6]. The largest part of the discrepancy results from the 
different system sizes: going from Cu+Cu to Au+Au at the same impact parameter 
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cuts the viscous reduction of v-i almost in half, to 39%. Viscous effects are much more 
dramatic for small collision fireballs than for larger ones, due to larger transverse 
velocity gradients. Accounting for different EOS used in [U [S] and [7J removes most 
of the remaining difference: replacing SM-EOS Q (with a first order phase transition) 
by EOS L (with a smooth cross-over fitted to lattice QCD data) reduces the viscous 
suppression of vi by about another third, to 28%. A final 5-10% relative reduction 
of the suppression results from replacing the " simplified I-S" treatment of [H [5] with 
the " full I-S" formalism used in [Jj . The combined effect of these changes is a viscous 
suppression of V2 of 25% in b — 7fm Au+Au at pr — 2 GeV/c, which agrees with [7J. 

6. Multiplicity scaling of elliptc flow and viscous entropy production 

Voloshin and collaborators pointed out (see [12] and references therein) that the 
eccentricity-scaled elliptic flow vije x scales with the final hadron multiplicity per 
unit overlap area, (l/S)(dN c h/dy). All dependence on system size, beam energy, 
and impact parameter can apparently be subsumed in this variable. The empirical 
scaling function rises monotonically with multiplicity density, approaching the ideal 
fluid prediction at high multiplicities but remaining well below it at low multiplicities. 

We investigated this "multiplicity scaling" for ideal and viscous hydrodynamics, 
for both vije x and the viscous entropy production fraction AS /So [B]. Even in the 
ideal fluid case we found weak scaling violations that can be traced back to the freeze- 
out condition which introduces an external scale that affects the fireball lifetime (see 
Fig. [I]). Perfect multiplicity scaling of v<i/e x in ideal fluid dynamics [13] can only be 
achieved only if decoupling happens late enough that the momentum anisotropy can 
fully saturate. In [6] we found this to be the case for very high collision energies. 

Scale breaking is more pronounced in viscous hydrodynamics and increases with 
rj/s. Experimentally seen weak scale breaking effects appear to be qualitatively 
consistent with viscous hydrodynamic predictions [6] but require more precise data 
for confirmation. 

Approximate multiplicity scaling is also found for the viscous entropy production 
fraction. AS /So depends strongly on the hydrodynamic starting time To. For 
To = 0.6fm/c and EOS L we find entropy production fractions ranging from about 
10% at high multiplicity densities to about 12-13% at low multiplicities; again smaller 
systems or more peripheral collisions exhibit stronger viscous effects. 
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